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1. Introduction 

NASA Contract NAS5-25898 entitled "Satellite Ranging Data Analysis" 
provided funding for the Department of Astronomy of the University of Texas at 
Austin under the direction of Dr. Peter J. Shelus, Principal Investigator, 
during the time interval 21 November 1979 through 20 August 1981. 

The original proposal spoke to six basic tasks which were to be undertaken 
as part of the contract. Subsequent negotiations toward an acceptable 
Statement of Wbrk gave high priority to two of the original six proposed 
tasks, i.e., 

- Task 1: Apply LAGEOS-derived polar motion components (x,y) to LLR 
data ; 

- Task 4: Using multi-station LLR data alone, compute values of Earth 
rotation continents (x,y) and UT1. 

Somewhat lower priority was assigned to two other of the proposed tasks, i.e., 

- Task 5: Relate LLR and LAGEOS reference frames; 

- Task 6: Compare Earth rotation components as derived using LAGEOS 
and LLR data against each other and against other techniques. 

Finally, the remaining two proposed tasks were not supported; 

- Task 2; Using residuals of LAGEOS ranging data with respect to a 
well-defined model, compute values of Earth rotation components 
(x,y) and UT1; 

- Task 3: Provide long-term calibration bench-marks for LAGEOS-derived 
uri values. 

2. MERIT Campaign Back-Up Analysis Center 

Under the subject contract, the university of Texas at Austin operated 
satisfactorily as a back-up Analysis Center (in addition to its role, under 
other funding, as a Computational Center) for the lunar laser ranging data 
type during the short MERIT Campaign which was held 1 August 1980 through 31 
October 1981. The term "back-up" refers to the circumstance that we were not 
obligated to respond in near-real-time. our responsibilities as back-up 

Analysis Center were to process LLR data gathered during the MERIT interval, 
subject it to suitable physical and mathematical modelling to produce a 
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UT1-UTC data product and to forward same to the MERIT Coordinating Center at 
the Bureau Internationale de l'Heure in Paris. Using a continually evolving 
system, designed for the extraction of Earth rotation information from the LLR 
data type (partially funded under this contract) , this task was satisfactorily 
handled in a routine manner. Hie data file transmitted to the BIH and the 
text report which accompanied that data transfer is included as Appendix A of 
this report. An informal comparison was made among the various LLR UT-1 data 
products at IAU Colloquium No. 63 at Grasse, France last spring. Hie 
comparisons were startling in their agreements. A formal report of this 
comparison will be included in the Proceedings of the Colloquium which should 
be available soon. Also, the McDonald Observatory LLR UT-1 MERIT data set is 
available in machine readable form either from the Principal Investigator or 
from the BIH. Since the results outlined above can employ polar motion 
parameters other than those supplied by the BIH in their monthly Circular D 
mailings, they satisfy, in part, contract requirements under Task 1 of our 
original proposal. 

Even though there was LLR data obtained from the Orroral Valley station in 

Australia, the data set forwarded to the BIH for the short MERIT Campaign 

contained only the results of analysis of McDonald Observatory data. Because 
of start-up problems at the Australian station during the interval of MERIT, 
this data set is still in the process of final filtering and data quality 
analysis. Part of this work of iteration to a final data set was performed at 
the Ihiversity of Texas under this contract in cooperation with Dr. Peter 
J. Morgan of the Australian station. Hie soft-ware system used and modified 
under this contract can process multi-station data and the work performed 

within the contract was included in our original statement of work, i.e., 

"...Using multi-station LLR data alone, compute values of Earth rotation 
components (x,y) and UT1..." and therefore satisfy, in part, contract 
requirements under Task 4 of our original proposal. 
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3. Intercomparison of IXR Results with Other Techniques 

?art of our efforts under this contract stressed some of the 
quick-turnaround and service oriented requirements of Earth rotation analysis. 
Therefore we strived to structure our algorithms and data processing systems 
to satisfy as many of these requirements as poiiible. Even though, by the 
expiration date of the contract, we were not able to transmit Earth rotation 
data products of the best possible quality because of temporary modelling 
deficiencies at the University of Texas, we were able to put into place all of 
the operational procedures for accomplishing same. These modelling problems 
concern mainly the lunar libration situation since we are totally dependent on 
external sources for this data product. Although efforts are underway both at 
the Massachusetts Institute of Technology and the Jet Propulsion Laboratory to 
supply us adequate modelling, such is not yet available. In spite of this, 
however, as already has been mentioned, preliminary intercomparisons of irr-1 
results have been made among the various analysis groups and observing 
techniques with very good results. Also, within days of contract termination, 
a file of some 18 months of ifTl-UTC information derived from McDonald 
Observatory LLR data was transferred by the University of Texas at Austin via 
the General Electric Mark III System to researchers at the U. S. Naval 

Observatory and the U. S. National Geodetic Survey for intercomparison with 
other research groups and various observing techniques. The data file 
transmitted to the USNO and the NGS is included as Appendix B of this report. 
The file can be supplied in machine-readable form to any and all persons who 
request it. 

Again, even though the data transferred to USNO and NGS consisted of 

information derived solely from McDonald Observatory LLR data, the entire 

system as designed and implemented can perform with multi-station data as 
well. Also, even though funding under this contract has expired, we remain in 
close contact, and are coordinating in whatever ways we can, with researchers 
from USNO and NGS to participate in the intercomparison among the various 
techniques in a scientifically meaningful way. Reports and publications 
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concerning these studie' are expected within the near future. The results 
reported in this section satisfy, in part, contract requirements under Task 6 
of our original proposal. 

4. Simultaneous LLR and LAGEOS Data Reduction 

Perhaps the most important output product which is a result of this 
contract and is also the task which received the greatest amount of effort 
during the term of this contract concerns the creation of a system which has 
the capabilities of simultaneously reducing LLR and LAGEOS laser ranging 
observations. Uhl ike most other studies of this type, we are not primarily 
interested in combining the Earth rotation results produced from the various 
techniques. Instead, we are interested in obtaining the Earth rotation 
parameters from the simultaneous reduction of the actual observational data 
itself. We feel that this is especially important when we consider both the 
lunar and the LAGEOS data types. Our study is trying to take advantage of the 
strong points of each data type and to eliminate, as much as possible, their 
weak points. We also stress the convenience for performing this task provided 
by the close proximity of the LLR analysis system within the UT Astronomy 
Department and the LAGEOS analysis system within the UT Department of 
Aerospace Engineering and Engineering Mechanics. 

As one might expect, the realization of such a reduction system has been 
fraught with many pitfalls. We have checked out this system with a rather 
large number of LLR data sets. As of the time of this final report, only two 
LAGEOS data sets have been made available to us for processing. However, 
based on only a modest checkout with respect to the LAGEOS data type, I feel 
that the system can accomplish all that had been planned and that the system 
can form the nucleus of a more complete system. It is hoped that additional 
funding will be provided by other means to make full use of this system. 

Regular status reports of the evolution of this system were made at the 
regular semi-annual meetings of the LAGEOS Investigators Meetings which were 
held at Goddard Space Flight Center. The definitive description of the 
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present state of this system was presented at the IAU Colloquiun No. 63 
entitled "High Precision Earth Rotation and Earth-Moon Dynamics, Lunar 
Distances and Related Observations" which was held in Grasse, France during 
22-27 May 1981. A preprint of the presented paper entitled "Earth Rotation 
from a Simultaneous Reduction of LLR and LAGEQS Laser Ranging Data" by Shelus, 
Zarate and Eanes is included as Appendix C of this report. 

The above-mentioned preprint covers the results of the efforts undertaken 
herein. The results presented in this section satisfy, in pert, contract 
requirements under Task 5 of our original contract. 

5. Administrative Summary 

As required under the terms of this contract, monthly financial summary 
reports, (NASA FORM 533M) were submitted to parsonell at GSFC (Code 269 and 
Code 942). Manpower expended during the tenure of this contract includes: Dr. 
Peter J. Shelus, Principal Investigator (3.6 mm); Dr. J. Derrai Mulholland, 
Research Scientist (2.2 mm); Wen-Jing Jin, Research Associate (3.0 mm); Gary 
Kern, Research Associate (0.9 mm); Alice Herzog, Clerk-typist (0.75 mm) . 
Trips were made by the Principal investigator to the semi-annual LAGEOS 
Investigators Meetings at GSFC and travel support to attend the IAU Colloquium 
Np. 63 was supplied by this contract. The foreign trip report submitted for 
that trip is attatched as Appendix D of this report. Finally, as of this 
writing, one formal publication has been submitted to the open literature, 
i.e., the paper by Shelus, Zarate and Eanes which is to appear in the 
proceedings of IAU Colloquiun No. 63. 

The principal Investigator wishes to that his colleagues at the University 
of Texas at Austin and elsewhere who collaborated on the efforts of this 
contract and who are too numerous to mention individually. He also wishes to 
acknowledge the support furnished to him by the GSFC Contracting Officer, 
Lauria A. Caria and the GSFC Technical Officer, C. C. Stephanides. 
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MERIT UT Results using McDonald LLR Observations 
at the University of Texas at Austin 

Peter J. Shelus and Nelson R. Zarate 

McDonald Observatory and Department of Astronomy 
University of Texas at Austin 
Austin, Texas 78712 USA 

These are UT resists computed from McDonald Observatory lunar laser 
ranging observations az the University of Texas at Austin. Since only 
single station data is being used, the analysis is similar to that per- 
formed by Shelus et al (1) . The lunar data set consists of some 63 normal 
points which represent some 600-700 Individual shot-by-shot ranges. Range 
residuals and partial derivatives are supplied by the standard LLR reduc- 
tion packages In regular use at the University of Texas. 

The LLR residuals are "post-fit", linearized residual and were 
obtained after a global parameter Improvement on 17 months of data approx- 
imately centered on the MERIT data set. Parameters in the global solution 
included linear, annual, and lunar nodal period terms in UT. Linearly 
interpolated values of BIH Circular D x, y, and UT1-UTC were used. 
Corrections to the Woolard nutation series, diurnal nutation terms in 
x and y, and diurnal tidal terms in UT as summarized by Williams (2) here 
applied. We are presently upgrading our LLR reduction system to more 
recent treatments of these corrections and the new IAU system of fundamen- 
tal constants. 

An additional note is in order concerning our results. Although a 
"window width" for observation selection is chosen similar to most other 
investigations (for instance, when computing 2-day averages, a window 
width of two days is chosen: all observations falling in that window enter 
into that particular solution), we do not necessarily slide the window 
one full window width before performing the next solution. Note that each 
individual solution is therefore not independent from neighboring 
solutions. In these results a window width of two days was used and we 
have slid this window by 0.04 days after each solution. A solution is 
performed only if there are at least three observations in a particular 
data set. After sliding the window, a new solution is performed only 
if the observation set has changed. 

(1) Shelus, P. J., Evans, S. W. , and Mulholland, J. D.: 1975, 

in Scientific Applications of Lunar Laser Ranging, ed. J. D. 

Mulholland, D. Reid el, Co, Dordrecht, P. 201. 

(2) Williams, J. G. : 1974, JPL Engineering Memorandum 391-592. 





obs 

ID 

726 
726 
726 
726 
726 
726 
726 
726 
726 
726 
72 6 
726 
726 
726 
72 * 
726 
726 
726 
726 
726 
726 
726 
72 6 
726 
72 6 
726 
726 
726 
72 6 
726 
726 
726 
726 
726 
72 6 
726 
726 
726 
72 6 
726 
72 6 
726 
726 
72 6 
72 6 
726 


A-2 




MODIFIED JD 
(.001 DAYS) 


UT I -070 

(.onnoisEC) 


STD. DEV. 

( .OOnOlSEC) 


44452697 

15767 

3 3 

44462860 

15737 

28 

44452987 

15722 

22 

44453115 

15671 

17 

44453330 

1562.0 

13 

44453502 

15586 

13 

44453672 

15552 

12 

44453389 

15522 

10 

4454016 

15493 

13 

44454184 

15462 

4 

44480431 

10836 

24 

44490683 

10773 

22 

44480807 

10760 

47 

44482233 

10286 

96 

44484951 

9756 

75 

44485132 

9705 

60 

44485313 

9690 

25 

44498762 

6684 

20 0 

44499 r ' 7 8 

6628 

65 

4 4 49 c : V/0 

6 56 5 

48 

44499698 

6506 

47 

44499913 

6493 

29 

44500251 

6420 

27 

44500473 

6348 

2 5 

44500593 

6317 

20 

44500976 

628 0 

18 

44501093 

6242 

20 

44501487 

6155 

27 

44501613 

6141 

22 

44501738 

6130 

25 

44502134 

6057 

1.8 

44502308 

6015 

30 

44502530 

5965 

40 

44502656 

5957 

33 

44502783 

5953 

39 

44502943 

5915 

48 

44503176 

58 55 

77 

44512537 

3188 

27 

44512758 

3135 

15 

44512921 

3097 

12 

44513046 

3106 

34 

44513172 

3086 

4 5 

44513385 

3009 

35 

44513555 

2971 

36 

44513783 

294 3 

42 

44513899 

2906 

65 


A-3 




725 

44514202 

2821 

41 

12 5 

44514545 

2779 

29 

72 5 

4451481 2 

2712 

33 

72 5 

44515193 

2545 

49 

725 

44515323 

2597 

39 

72 5 

44515514 

2551 

42 

72 5 

4451 5845 

2500 

25 

72 5 

4451595’ 7 

2481 

41 

725 

44534242. 

-2057 

42 

72 5 

44534457 

-2133 

25 

72 5 

44534527 

-2187 

27 

725 

44534753 

-2210 

30 

725 

44534925 

-2252 

38 

72 5 

448V734P 

-3112 

15 

725 

44538044 

-3300 

22 

72 5 

44538738 

-3490 

40 


APPEtDIX B 



DBS 

ID 

72* 

72 * 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72 * 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 

72* 


B-l 


ORIGINAL PAGE 13 
OF POOR QUALITY 


MODIFIED JD 
(.001 DAYS) 

4^302**0 

44302839 

4430297* 

4430.31*3 

44303303 

44303*32 

44303792 

44303919 

4430404* 

4430420* 

44304429 

44333031 

4433479* 

4433*487 

44335720 

44335880 

4433*008 

4433*180 

443**818 

4435*954 

44357241 

443*4475 

443841*2 

44384391 

44384548 

44384*74 

44384800 

4438501* 

44385188 

44385413 

44385525 

44393120 

44393233 

443934** 

44393*9* 

44393812 

44395205 

44395900 

4439*309 

4439*553 

4441*5*2 

4441*740 

44417.739 

44417452 

44417729 

44417925 

44418257 

44418583 


UT1-UTC 
( .00001SEC) 

4872* 
49757 
48707 
48*28 
48527 
48477 
48439 
48405 
48340 
48301 
48245 
40871 
4 0522 
40234 
4 0200 
4015* 

4 01*4 
40127 
34751 
34702 
34593 
32829 
28079 
28047 
28012 
279*9 
27957 
278** 

278 28 
277*0 
27*95 
2*308 
2*29* 
2*2*7 
2*215 
2*227 
25953 
25*85 
75517 
2547* 
21*09 
2158 2 
21519 
21480 
21407 
21400 
21388 
21270 


STD. DEV. 

( . 00001SEC) 

1 . 1 * 

110 

55 

54 
51 
45 

37 
27 

4 

5 
8 

13 

10 * 

2 * 

33 

39 
3 3 
43 
94 

38 
2 

21 
17 
27 
2 2 
7.0 
21 
3 4 

40 
5* 

25 

50 

34 
29 
42 
27 
21 
25 
88 

123 

1 

1 

13 

33 

55 
94 
77 
84 



72.6 

4441R779 

7 2 6 

44418975 

726 

44419297 

726 

44419628 

726 

44419818 

"7 2 6 

44420007 

726 

4444078? 

72.6 

44443789 

726 

44444112 

726 

44444607 

726 

44444727 

726 

44444894 

726 

44445131 

726 

44449364 

726 

44449595 

72 6 

44449714 

726 

44450111 

726 

44460361 

726 

44452041 

72 6 

44452177 

726 

44452468 

72.6 

44452697 

726 

44452860 

726 

4 4452987 

726 

44453115 

726 

44453330 

726 

44453502 

726 

44453672 

72 6 

44457RR9 

726 

44454016 

72 6 

44454184 

72 6 

44480421 

72 6 

44480683 

7 2.6 

44480807 

726 

44463233 

7 7.6 

44484051 

72 6 

44485132 

726 

44485313 

726 

44498762 

72 6 

44499078 

77.6 

44499580 

726 

44499698 

726 

44499913 

726 

44500251 

72 6 

44500473 

726 

44500593 

726 

44500926 

726 

44501093 

77.6 

*4501487 

72 6 

44501613 

726 

44501738 


2 

of'SL Nal 

0F POOR 

21200 

52 

21214 

73 

21222 

57 

21145 

65 

2108 s 

57 

21087 

54 

17596 

109 

17153 

100 

17202 

100 

17075 

95 

17032 

69 

17037 

78 

16983 

136 

.16405 

11 

16354 

13 

16327 

20 

16252 

.12 

16207 

2 

15997 

56 

15921 

29 

15839 

23 

15767 

33 

15737 

28 

15722 

22 

15671 

17 

15620 

13 

15586 

13 

15552 

12 

15522 

10 

15493 

13 

15462 

4 

10876 

2 4 

10773 

22 

10760 

47 

10286 

96 

9756 

75 

9705 

60 

9690 

25 

6684 

200 

6628 

65 

6 56 5 

48 

6506 

47 

6493 

29 

6420 

27 

6348 

25 

6317 

20 

628 0 

18 

624 2 

20 

6155 

27 

6141 

2 2 

6130 

2 5 



B-3 


725 

44502134 

5057 

725 

44502308 

8015 

725 

44502530 

59 5 5 

725 

44502555 

5957 

725 

44502783 

59 53 

725 

44502943 

5915 

72 5 

44503175 

58 55 

725 

44512532 

3188 

725 

44512758 

3135 

725 

44512921 

3097 

725 

44513045 

3105 

725 

44513172 

3085 

72 5 

445133R5 

3009 

725 

44513555 

2971 

725 

44513783 

2943 

725 

44513899 

2905 

725 

44514302 

2321 

725 

44514545 

2779 

72 5 

44514813 

2712 

725 

44515193 

2545 

72 5 

44515323 

2597 

725 

44515514 

2551 

725 

44515845 

2 50 0 

725 

44515957 

2 481 

72 5 

44534242 

-2057 

”725 

44534457 

-2133 

725 

44534527 

-2187 

725 

44534753 

-2210 

725 

44534925 

-22 57 

72 5 

44537348 

-3112 

72 5 

44538044 

-3300 

725 

4453873R 

-3490 

72 5 

44500139 

-18525 

725 

44500857 

-18755 

72 5 

44832553 

-25353 

725 

44548708 

-30112 

725 

44548845 

-30085 

725 

44549135 

-30132 

725 

44549354 

-30.159 

72 5 

44549525 

-30207 

725 

44549555 

-30219 

725 

44549783 

-70225 

725 

44549945 

-30253 

725 

44550174 

-30305 

725 

44550853 

-30383 

72 5 

44551551 

-30530 

725 

44552242 

-30557 

72 5 

44552473 

-30587 

72 5 

44552535 

-30731 

725 

4 45 52. ^ 5 5 

-30593 

72 5 

44852940 

-30733 


18 

30 

40 

33 
39 

48 
77 
37 
15 
12 

34 
45 

3 5 

35 

4 2 
55 

41 

29 

3 3 

49 

39 

42 
25 
41 

4 2 
25 
27 

30 
3R 
15 
22 

40 

50 

A 


59 
43 
29 
31 
2 5 
21 
20 
15 
17 
24 
14 
8 3 
95 
57 
7 2 
55 
92 


B-4 


ORIGINAL PAGE IS 
OF. POOR QUALITY. 


726 

4465977R 

-32276 

12 

726 

44659960 

-32313 

29 

726 

44660095 

-32315 

33 

726 

44660292 

-32356 

50 

726 

44660510 

-32368 

4 

72 6 

44677449 

-37303 

51 

726 

44677622 

-37344 

36 

726 

44677793 

-37387 

14 

72 6 

44680545 

-38093 

105 

726 

44681233 

-38167 

42 

726 

44681923 

-38342 

64 

726 

44682040 

-38387 

4 2 

72 6 

44682272 

-38438 

31 

72 6 

44682497 

-38527 

29 

726 

44682617 

-38532 

20 

72 6 

44686381 

-39501 

2 

726 

44705483 

-44548 

36 

726 

44719481 

-48285 

26 

72 6 

44719712 

-48328 

21 

726 

44719876 

-48362 

20 

726 

44720003 

-48392 

30 

726 

44720294 

-4' 4 61 

21 

726 

44720497 

-48512 

22 

72 6 

44720798 

-48539 

20 

726 

44736797 

-52572 

11 

72 6 

44737482 

-52733 

1 '6 

726 

44738206 

-52695 

3 6 

726 

4473^3.49 

-53110 

18 

726 

44746672 

-54981 

8 2 

72 6 

44746866 

-55128 

106 

726 

44747003 

-55088 

63 

726 

44747182 

-55133 

90 

726 

447474^2 

-55183 

121 

726 

44748130 

-55292 

102 

72 6 

44748820 

-55426 

79 

726 

44749498 

-55484 

21 

726 

44749738 

-55585 

67 

72 6 

44749853 

-55688 

28 

726 

44750261 

-55716 

71 

72 6 

44750509 

-55776 

104 

726 

44751176 

-55730 

162 

726 

44751299 

-55933 

98 

726 

44751580 

-55968 

145 

726 

44781833 

-55945 

388 

726 

44751949 

-55602 

732 

726 

44766470 

-59352 

19 

726 

44766660 

-59387 

14 

726 

44766943 

-59471 

2 2 

72 6 

44767149 

-59513 

20 

726 

44767351 

-59552 

28 

726 

44767629 

-59617 

4 5 



B-5 


page is 

"** QUALrfi 


726 

44767760 

-59656 

30 

72 ^ 

4476703 ft 

-59704 

41 

726 

44768173 

~ 59 7 73 

44 

7 ?.6 

44768044 

-60003 

on 

726 

44777465 

- 612.81 

46 

72 6 

44706771 

35777 

17 

7 26 

44796803 

35758 

7 

726 

44797126 

3 5713 

19 

726 

4479737 ft 

35694 

23 

726 

44797538 

35684 

22 

726 

44797663 

35688 

19 

726 

44797785 

35664 

18 

726 

44707946 

- .35649 

17 

72 6 

44798197 

35626 

22 

726 

44798431 

35628 

21 

72 6 

44798606 

3 550 5 

28 

726 

44798744 

35553 

53 

72 6 

44798932 

35520 

78 

72 « 

44799261 

3 5613 

142 

726 

44790585 

35539 

24 ? 

72 6 

44806721 

34873 

149 

726 

44806915 

34862 

87 

72 6 

44807058 

34838 

45 

72 6 

44807241 

3 4 7 6 2 

27 

726 

44807472 . 

34751 

2 5 

726 

44807725 

34719 

27 

726 

44807907 

34719 

32 

726 

44808040 

34772 

38 

726 

44808220 

34695 

54 

726 

448.37006 

30038 

144 

726 

44837235 

29912 

74 

726 

44837528 

29807 

47 

72 6 

44837755 

2974 5 

46 

72 6 

44837871 

29692 

9 


APPENDIX C 



C-l 


Earth rotation from a simultaneous reduction of 
LLR and LAGEOS laser ranging data 


Peter J. Shelus , Nelson R. Zarate and Richard J. Eanes + 
^McDonald Observatory and Department of Astronomy 
+Department of Aerospace Engineering and Engineering Mechanics 
University of Texas at Austin 
Austin, Texas 78712 USA 


I. INTRODUCTION 

As the techniques of lunar and artificial satellite laser ranging 
mature, emphasis is being placed upon the use of these observations 
to monitor the Earth's rotation. It is important to note, however, 
that at the present time neither technique alone can furnish al l 
three components of this rotation to an accuracy which surpasses 
those results obtained from classical techniques. In the case of 
LAGEOS laser ranging, unmodeled secular orbital effects couple with 
axial Earth rotation in such a way that these effects are not 
separable in the analysis of those observations. In the case of 
lunar laser ranging, observations have been regularly available only 
from a single station for the past ten years or so with the result 
that a change in latitude along the McDonald Observatory meridian is 
not separable into "he ordinary (x,y) components of polar motion. 
The main purpose of this paper is to present the first stages of an 
investigation to combine LAGEOS and lunar laser ranging 
observations. It is hoped that the proper implementation of such a 
process might eliminate the shortcomings inherent in each technique, 
while accentuating the advantages of each. This has the potential 
of producing all three components of the Earth's . rotation to an 
accuracy and precision which is compatible with the present 
observational uncertainties. 

II. DATA AND MODEL COMPATIBILITIES 

As is the case in all investigations which seek to combine two 
different observation types, a great deal of ground work must first 
be laid before the data synthesis can be begun. Care must be taken 
to insure that the various data types to be combined are totally 
compatible and consistent with one another. Not only must such 
mundane matters as units and formats be unambiguously defined, but 
standards for such things as reference frames, theoretical and 
empirical models as well as fundamental and derived constants nust 
be strictly and totally adhered to. 



At the University of Texas at Austin, two parallel efforts are 
underway to obtain Earth rotation information from artificial 
satellite and lunar laser ranging observations. That within the 

Department of Aerospace Engineering and Engineering Mechanics has 
been concerned with the artificial satellite analysis; that within 
the Department of Astronomy and McDonald Observatory has been 
concerned with the lunar analysis. Although each has been 
performing its tasks completely independently of the other, each 
uses the extensive computing facilities of the University's main 
computer systems. This happy circumstance has alleviated many of 
the problems associated with the transfer of data and information 
between independent reduction and analysis systems. Also, the close 
proximity of the personnel of both groups assures precise 
communications and thereby has eased the reference model 
compatibility problems. 

As might be expected, our initial efforts have been applied to 
testing algor ithms and applying them to the LAGEOS and LLR data sets 
which were obtained during the short MERIT campaign which ran from 
August through October of 1980. The lunar data set consists of some 
63 normal points which represents some 600-700 individual lunar 
laser ranging observations. Specific information about thi^ data 
set can be found in the MERIT Campaign Report which should be 
generally available from the Bureau Internationale de l'Heure in 
Paris. The LAGEOS data set contains in excess of 20,000 individual 
LAGEOS ranges and will be described elsewhere. In both cases, our 
analysis efforts are concerned with range residuals and partial 
derivatives which are supplied by the standard LLR and LAGEOS 
reduction packages which have been in regular use over the past few 

years at the University of Texas. Although, to the best of our 

knowledge, the current data sets are internally consistent and 
compatible, additional checks will be made continuously throughout 
the course of the total investigation to preserve and/or extend this 
integrity. 

III. ANALYSIS TECHNIQUES 

Because of the short term nature of the effects being sought by this 
investigation, our "observational equation" is a simple one and, for 
the LLR case, is similar to that presented by Stolz and Larden 
(1976), i.e., 

pQ-fc = r [ sin cos d cos(X-H) - cos \ sin d cos $ ] x 

-r ( sin <}> cos 6 sin(X-H) - sin X sin 6 cos § ] y 

+r cos <j> cos d sin H J(UTl-UTC) 

where Pq-'Pc * s t ^ ie ran 9 e residual (observed minus computed) ; r is 
the radius of the Earth and X is its east longitude; K is the local 
hour angle of the retroref lector and d is its declination; x, y, and 
ci(UTl-UTC) are improvements to the nominal values of these Earth 
rotation parameters. Although the above expression is that which is 
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specifically used for LLR, a similar one has been used for the 
LAGEOS case « 

This investigation is seeking values for the Earth rotation 
parameters averaged over 5 day intervals or less. In the case of 

LLR these short-term effects are well-separated from any unmodelled 
long term effects because it is believed that all short term (less 
than two weeks or so) lunar orbital and librational effects down to 
the few centimeter level are known. This is, of course, not yet the 
case for LAGEOS and it is certainly recognized that analysis efforts 
to extract orbital information from the LAGEOS data also extract 
axial Earth rotation information from that data, thus decreasing 
one's ability for obtaining accurate UT1-UTC information from this 
data type. It is believed that this study is the first attempt to 
obtain Earth rotation parameters by the simultaneous reduction of 
LLR and LAGEOS data at the observation level. 

The LLR residuals which were used in this study are "post-fit" , 
linearized residuals having been obtained after a normal global 
parameter improvement run on some 17 months of data centered 
approximately on the MERIT data set. Parameters in the global 
solution runs include linear, annual and lunar nodal period terms in 
UT. Linearly interpolated values of smoothed BIH Circular D x, y, 
and UT1-UTC which were modified by corrections given by Williams 
(1974) based on McClure (1973) have been used r Also used was the 
Woolard (1953, 1959) nutation series as modified by Melchoir (1971). 
Simple checks have shown slight differences with the Wahr (1980) 
nutation series and the Yoder et al ( 1 L>8 1 ) treatment of UT diurnal 
tidal terms, although we are presently upgrading our LLR reduction 
systems to these more recent treatments and the new IAU system of 
fundamental constants. 

The LAGEOS residuals were computed with the model used to generate 
the LAGEOS long-arc trajectory designated LLA80.il. The gravity 

field used was LGM80.11.1; this geopotential is a preliminary 
LAGEOS-derived adjustment to the GEM10 field. The model includes 
the Wahr nutation series, the short period variations in UT from 
Yoder et al (1981), and BIH Circular D smoothed values for polar 
motion and UT1. The station positions were the LAGEOS-derived set 
designated LSC00.11. Orbit initial conditions were estimated from a 
sampled set of LAGEOS observations from 16 sites over the 124 day 
period from 30 June to 31 October 1980 (MJD - 44420-44543) . The 
full set of data contained 508,000 observations while the sampled 
set, obtained by requiring that no two observations from any one 
site be less than one minute apart, contained 22,000 ranges. The 
unweighted RMS of the post-fit residuals was 0.42m. The estimated 
"single-shot" precision was 0.25m when averaged over all of the 
laser systems involved. 

The remaining unmodelled long-period variations in the LAGEOS 
orbital elements were removed by smoothing the element residuals 
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from LLA80.il with a Vondrak filter using £=1.0E-06 {half power at 
60 days) , Because of the high correlation of errors in UT and 
errors in the LAGEOS orbit node, this empirical adjustment to the 
LAGEOS orbital elements effectively filters a portion of any signal 
present in UT1-UTO. The small correlation of polar motion 

components x and y to the orbital elements implies that they are 
only slightly affected by the empirical adjustment. As the LAGEOS 
dynamical model matures the use of an empirically corrected orbit 
will be discontinued. 

IV. NUMERICAL RESULTS 

Using an observational equation of the type give in Section III, we 
have computed observational residuals and partial derivatives using 
standard lunar and LAGEOS data analysis packages. Several of the 
initial solution attempts are being reported here. To assess the 
solution algorithms of this package the first solution run was 
performed to obtain UT1-UTC and a constant bias from LLR 
observations alone. Since only single station LLR data is being 
used in this study the analysis is similar to that performed by 

Shelus et al (1976) . Ihe second solution run was performed to 

obtain x, y f and UTl-UTC estimates from LAGEOS observations alone. 
Only those LAGEOS observations which were close in time to LLR 

observations were used (a full analysis of the LAGEOS-only results 
is beyond the scope of this paper) . Each of these two runs provided 
results which were similar to those obtained from analyses performed 
independently of this study. The results, which give deviations to 
BTH Circular D 5-day smoothed values, can be seen for UT1-UTC in 
Figures 1 for the LLR-only case and for x, y, and UT1-UTC in Figure 
2 for the LAGEOS-only case. 

From an examination of these figures we see that our initial 

expectation that UT1-UTC "power” has been lost from the LAGEOS 
observations is confirmed since the deviations from BIH values for 
UTl-UTC are much smaller from the LAGEOS-only results than from the 
LLR-only results. This assumes, of course, that the LLR-only 
results are "correct". Having confirmed our expectations, we next 
proceed to the next step whereby we may "tie” ’ the short-term 
signature from the LAGEOS data type to the long term signature from 
the LLR data type by attempting simultaneous solutions. Figure 3 
shows the results for our first such attempt. In this case we have 
opted to only consider the x and y partial derivatives (not UT1-UTC) 
from the LAGEOS data set simultaneously with all three partial 
derivatives (x, y, and UT1-UTC) from the LLR set. All observations 
going into the solutions are given equal weight in spite of the 
overwhelming amount of LAGEOS data with respect to the LLR data. 
Ihe signature for the UT1-UTC results are similar to the LLR-only 
results, and the signatures for the x and y results are similar to 
the LAGEOS-only results, as would be expected. 

A very important sidelight of this investigation surfaces from our 



processing of the observations in a manner different from most. 
A1 chough a “window width” is selected similar to most other 
investigations (for instance, when one is computing tw day 
averages, one chooses a window width of two days; all observations 
which fall through that window are allowed to enter into that 
solution) , we do not necessarily move our window one full window 
width before performing the next solution. We feel that this 
technique can give a more complete representation of the information 
content of each data set. However it does have the drawback that 
each individual solution run is not completely independent from 
neighboring solutions. In all of the results winch are presented 
here we have used a "window width" of two days and have "slid" this 
window by 0.04 days after each solution. A solution is performed 
only if there are at least three LLR observations in a particular 
data set. After sliding the window, a new solution is performed 
only if the LLR data set has changed. 

The results presented in Figure 3 are not satisfying from several 
points of view. First, there were no attempts made to normalize the 
effects of each observation type through proper weighting 
parameters. As has already been mentioned, the LLR data are normal 
points while the LAGEOS data are shot-by-shot data. A far more 
serious objection arises from the fact that only a very weak tie is 
established between the tv» data types because the UT information 
from the LAGEOS data has been ignored and only single station LLR 
data exists. A crude attempt at normalization was made for the 
fourth solution run (Figure 4) wherein the third solution was 
performed aga:n except that the LLR normal points were given a 
weight 5.0 with respect to unity for the LAGEOS shot-by-shot points. 
As might be expected, Figures 4 and 5 are quite similar. 

V. DISCUSSION 

Although the results presented here are preliminary, they are 
indicative of the great progress which has been realized recently at 
the observation by observation level in the combination of LAGEOS 
and LLR results for Earth rotation. Each technique is certainly 
mature enough that consistency and compatibilty between such 
different data types has been accomplished. The presence of such a 
tw-pronged analysis effort opens the door to a more proper and 
satisfying data synthesis. Our next steps will progress to more 
realistic ties between the two data types. This will entail using 
the LLR results to help separate the unmodelled orbital effects of 
LAGEOS from axial Earth rotation instead of merely ignoring the 
effects of UT1-UTC in the LAGEOS data. Simultaneously the x and y 
results from LAGEOS will be used to improve the LLR results. Chce 
separation is obtained the short term LAGEOS results will be 
"anchored" by the long-term LLR results, thereby giving the UT1-UTC 
parameter the same significance and resolution as the x and y 
parameters. 
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Further progress will be also accomplished by a farther 
investigation, of the relative weighting schemes for LAGEQS versus 
LLR data to more reasonably combine normal point and shot-by-shot 
data. It may be also attempted to work with UX\ shot-by-shot data 
and/or LAGEOS normal point data to obtain this next level of 
compatibility. 
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Figure Captions 



Figure 1. Differences in UTl-UTC determined by LLR with respect 
to Bill Circular D 5-day smoothed values. 


Figure 2. Differences in x, y, and UTl-UTC determined by LAGEOS 
with respect to BIH Circular D 5-day smoothed values. 


Figure 3. Differences in x, y, and UT1-UTC determined by LLR and 
LAGEOS with respect to BIH Circular D 5-day smoothed 
values (LAGEOS sampled shot-by-shot data and LLR normal 
point data equally weighted). 


Figure 4. Differences in x,y, and FTl-UTC determined by LLR and 
LAGEOS with respect to BIH Circular D 5-day smoothed 
values (LLR normal point data weighted by a factor 5 
with respect to LAGEOS shot-by-shot data). 
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THE UNIVERSITY OF TEXAS AT AUSTIN 

COLLEGE OF NATURAL SCIENCES 
Department of Astronomy 

AUSTIN, TEXAS 78712 D?te: 2 June 1981 

MEMORANDUM 

To: 

From: 

Subject: 


C. C. Stephanides (Code 904) and L. E. Walker (Code 269) 
Peter J. Shelus K 

Foreign Travel Trip Report (Contract NAS5-25898) 


During May 18-21 I attended the Project Merit Workshop which 
was held in Grasse, France. This workshop was held under the 
auspices of the 1AU/IUGG Joint Working Group on the Rotation of the 
Earth. 

Informal sessions were held among participants on Monday 
morning. On Monday afternoon brief formal summary reports were 
presented by the directors of the various observation, analysis, 
computing and communication centers. These summary reports were 
mainly concerned with what went right and what went wrong during 
the short MERIT campaign. A short discussion on changes which 
might be recommended to insure a successful long campaign followed. 

The entire day Tuesday was dedicated to expanded reports on 
the short MERIT Campaign from the various technique coordinators 
and their sub-coordinatcrs . During these sessions I presented a 
report on the McDonald Observatory participation in the short MERIT 
campaign . The report summarized the data acquisition and data 
communications tasks which were accomplished. 

Formal reports continued into Wednesday. I also attended 
meetings of the MERIT Standards Committee chaired by W. G. Melbourne, 
the lunar laser ranging sub-group for MERIT chaired by J. D. Mulholland 
and the EROLD Steering Committee also chaired by Mulholland. Topics 
related to LLR were related to the modelling standards which would 
be used during the regular MERIT campaign. Most crucial to LLR 
are of course the lunar orbit and rotation models. Along those lines 
I was assigned the action item of communicating with appropriate 
individuals and submitting a final report to the MERIT Standards 
Committee by October 1981 with LLR recommendations. 

A General discussion on the short MERIT campaign and a review 
of aims and tasks for the regular campaign occupied us on Thursda '. 

Of most importance for LLR. was the decision to separate LLR and SLR. 
During the short campaign Aardoom and Mulholland were subcoordinators 
for the artificial satellite and lunar techniques respectively under 
Silverberg. For the main campaign 0. Calame (of CERGA) will be the 
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coordinator for the LLR technique while J. La t timer (of SAO) will 
be the coordinator for the SLR technique. There will be no "overall" 
coordinator for the laser techniques. A final decision was to 
hold the regular MERIT campaign from September 1, 1982 through 
October 31, 1983. 

In summary, it was generally accepted that most of the goals 
of the short MERIT campaign were successfully attained. Each of the 
techniques were able to observe, reduce and transmit Earth rotation 
information to the Coordinating center at BIH in a timely manner. 

The use of the electronic transfer of information was of great 
importance. However, a significant amount of cooperation and 
communication will need to be exercised to assure the efficient 
and economical use of the system during the up-coming long MERIT 
campaign. A sub-group consisting of Peter Morgan, Jim Lattimer, 
Martine Feissel and myself were assigned the task of looking into 
these matters. 

Finally, it was acknowledged that the short MERIT campaign also 
formed the incentive for many of the techniques to come "up to speed" 
earlier than they might have under ordinary circumstances. The 
effect of the Chinese classical data did much to improve the 
results of that technique. We are all looking forward to an equal 
success of the regular campaign. 

On Friday, May 22 and Monday through Wednesday, May 25-27, I 
attended the sessions of IAU Colloquium No. 63 entitled "High Precision 
Rotation and Earth-Moon Dynamics, Lunar Distances and Related 
Observations" at Grasse, France. 

Friday served as a transition day as the MERIT Workshop wound 
down and the Colloquium commenced. A total of 12 papers were pre- 
sented on Friday which summarized the scientific results which were 
obtained by the various observational techniques for Earth Rotation 
monitoring during the short MERIT technique. I presented a paper 
as this session which was entitled "Earth Rotation from a Simul- 
taneous reduction of LLR/LAGEOS Data". 

On Monday 4 papers were presented in the morning which supplied 
historical and complimentary material related to the concepts of 
Earth rotation monitoring. The afternoon was spent in a meeting of 
the EROLD Steering Committee. This meeting was essentially a 
continuation of that begun by the LLR group during the MERIT Workshop 
the previous week. Mainly discussed were the relations, if any, 
which should be eliminated, maintained or inaugurated between the 
EROLD and the MERIT campaigns. 
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Tuesday and Wednesday saw a total of 18 papers being presented 
on many of the aspects of lunar dynamics. These papers ran the 
gamut from energy dissipation, tidal friction and relativistic 
perturbations to analytical and semi-analytical theories of the 
Moon as well as reference frame determinations. 

I returned to the U.S. on Thursday, 28 May 1981. 


PJS/ah 



